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Abstract 

We compute the particle currents induced on a bubble wall background at finite temperature 
in a model with CP violation in the Higgs sector. Using a field theory approach we show that 
fermionic currents arise at one loop, so that a suppression factor 0(ht<f)/TTT) 2 with respect to 
previous computations is found. The contributions to the Higgs currents are also derived and 
their relevancy for the spontaneous baryogenesis mechanism is discussed. 
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The possibility of generating the baryon asymmetry of the Universe during the electroweak 
phase transition has received much attention in the last years § . In the usual scenarios the 
transition is required to be first order and to proceed via nucleation of bubbles of the broken 
phase in the unbroken phase. The necessary departure from thermal equilibrium then can take 
place inside or in front of the walls of the expanding bubbles. 

Two different limits have been investigated in the literature. In the case of thin bubble walls, 
the asymmetric (in fermion numbers) reflection of particles off the bubble wall is the dominant 
effect. The induced fermion number flux is then reprocessed into a baryon asymmetry by the 
anomalous, (B + L)-violating, sphaleronic transitions in the unbroken phase 

In this letter we will focus on the opposite limit of thick bubble walls. Indeed, this is thought 
to be the relevant one if the phase transition is weakly of the first order, as it seems to be the 
case for the standard model 0, |5| and the minimal supersymmetric extension of it (MSSM) [FJ. 
If the only source of CP violation is the phase in the Cabibbo-Kobayashi- Maskawa matrix, it 
seems very hard to generate any baryon asymmetry in this limit. On the other hand, if CP 
violation is present in the Higgs sector (which requires at least two Higgs doublets), the so 
called spontaneous baryogenesis mechanism can be invoked 0. Both in the case of explicit 
@ and spontaneous || CP violation, a space-time dependent relative phase 5(x) between the 
two Higgs fields is turned on inside the bubble wall. In order to analyze the effect of this 
complex space-time dependent background on particle densities, a rotation on the fields can be 
performed to make the Yukawa couplings real. As a consequence, a derivative coupling of the 
form 

C int ~ d^J>*, (1) 

where J M is the current corresponding to the rotationF], is induced from the kinetic terms. 

In the original paper 0, only the time derivative of the phase 5 was taken into account. In 

this approximation, 5 acts as an effective chemical potential (usually called 'charge potential'), 

and the particle densities are perturbed to nonzero values rij ~ 5 T 2 , where is the charge 

of the i-th particle under the given rotation. In presence of such a charge potential, and of 

baryon number violation, the thermodynamical evolution of the system adiabatically leads to 

a non- vanishing baryon asymmetry. This approximation has been improved in ref. 0, where 

also the spatial derivatives of 5 in eq. (|l|) were taken into account, and the role of particle 

order to avoid anomalies, in the original paper fif] and in U the rotation was taken to be the hypercharge. 
In fact, this can only be done if the same higgs doublet couples both to up and down type quarks. In the models 
in which two different Higgs doublets are coupled to the up and down type quarks (as is the case for the MSSM) 
it is not possible to remove the relative phase S from all the Yukawa couplings by making a hypercharge rotation 
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diffusion was discussed. 

Nevertheless, we believe that using the interaction term in eq. ([!]) as a starting point to 
compute the perturbations to the thermal averages presents some problems. 

First of all, the proportionality between an individual particle density and its hypercharge 
is a direct consequence of the hypercharge rotation made to make the Yukawa couplings real. 
Making a different rotation, the proportionality would of course drop. For instance, one could 
rotate the right handed top and leave the left handed one untouched, absorbing the anomaly 
by a proper rotation of the light fermion fields. In this case, the perturbations to be put in the 
system of kinetic equations describing the system would be different, so that in principle the 
final value for the baryon asymmetry might come out to be dependent on the rotation that has 
been made. 

Secondly, since the phase 5 is communicated from the Higgs to the fermion sector through 
the Yukawa interactions, any perturbation in the fermion densities rii should vanish in the limit 
of zero Yukawa couplings hi. Also, they should vanish in the limit of zero vacuum expectation 
value for the Higgs fields Hi(x) = Vi(x) exp[z0j(x)] because no spontaneous CP violation is 
present in the Higgs sector in this limit. Naively, one could then expect a suppression factor 
of order (h 2 vf(x)/T 2 ), where hi is the relevant Yukawa coupling, for the perturbations in the 
fermionic particle number with respect to the original result. Since we are interested in regions 
of the bubble wall where sphalerons are still active, i.e for values of Vi(x)/T typically smaller 
than one, then the above mentioned suppression factor might be crucial. 

The ultimate reason why these suppressions do not appear in the original treatment, is that 
considering eq. (|l|) as the only effect of the background is equivalent to perturbing around the 
Higgs field configuration 8(x) = 0, Vi(x) 7^ 0, which is not a solution of the field equations. In 
other words, it is equivalent to disentangling 5(x) from Vi(x), whereas from the field equations 
one can see that 8^5 (x) vanishes as Vi(x) 2 for vanishing Vi(x). 

The purpose of this letter is to compute the averages on the bubble wall background, 
both for fermions and Higgses, making use of a field theoretical approach. In this way we are 
able to treat the background consistently and to recover the expected suppression factors. We 
wish to stress that what we call here n$ are just the perturbations induced by the CP violating 
background. They should be used as source terms for the departure from equilibrium in the 
equations describing dynamical processes, like gauge and Yukawa interactions, baryon number 
violation and particle diffusion. 

Our starting point is the finite temperature generating functional for the 1PI Green's func- 
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tions with insertion of an operator O(z) (in the following O(z) will represent a particle current) 

T A(x)} = W [Ji(x), A(x)} -Y,J cPxJjix^fa), (2) 

5 

where §\{x) are the classical fields of the theory and Ji(x) the corresponding sources, while 
A(x) is the source for the operator 0(x). 

The quantity we are interested in is the expectation value of the operator O(z) on the 
background given by the fields which we will specify later, 

1 5T A] 
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We can expand the functional O [^{x)] (z) in a power series of $^ 
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where the coefficients of the expansion are the n-point 1PI Green's functions with one insertion 
of the operator 0(z) computed in the unbroken phase 
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5 n+1 T [$?, A] 
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(5) 



The relevant background $f(a;) for us is given by the bubble wall, which is a solution of the 
field equations of motions 

5T [$?, A = 0] 







(6) 



with appropriate boundary conditions. 

In order to be more specific, we consider the extension of the standard model with two Higgs 
doublets where the up and down type quarks couple to different Higgs doublets. We consider 
here the case of spontaneous CP violation like that considered in ref. [Q for the MSSM. 

The most general tree level effective potential is given by 

V = m 1 2 \H 1 \ 2 + m 2 2 \H 2 \ 2 -(m 3 2 H 1 H 2 + h.c.) + X 1 \H 1 \ 4 + X 2 \H 2 \ 4 

+ A 3 |#i| 2 |# 2 | 2 + A 4 |#i# 2 | 2 + [A 5 (#i# 2 ) 2 + MH^HrHz + \ 7 \H 2 \ 2 H X H 2 + h.c] , (7) 

where all the couplings are now assumed to be real. Note that the potential V depends only 
on the phase 5 = 9 X + 9 2 , whereas the orthogonal combination represents the gauge phase. 

We assume that the parameters of the Lagrangian and the temperature are such that, when 
the loop corrections are considered, the potential assumes a double well shape and a expanding 
bubble solution exists for the complete equation of motion given by eq. (§). 
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Since the interesting dynamics for baryogenesis takes place in a region close to or inside the 
bubble wall, we approximate it with an infinite plane travelling at a constant speed v w along 
the z axis. 

Considering eq. (H) at the tree level and including the one loop corrections to the effective 
potential, one can see that a solution exists for which all the background fields are vanishing 
except for the neutral Higgses. Their asymptotic behaviour is such that Vi(z) — > for z — > — oo 
and Vi(z) — > vf for z — > +00, where vf are the finite temperature values for the Higgs fields in 
the broken phase. Moreover, due to the dependence on the phase 5 of the potential, the two 
Higgs phases 6*1,2 are space-time dependent, and their evolution is such that 

vl{z)& t e x (z)=vl(z)& i e 2 {z). (8) 

This solution is the starting point for our expansion in eq. (^), from which we will now 
see that a non zero contribution to the neutral Higgs currents already exists at the tree level, 
whereas the currents for fermions and charged Higgses appear as loop effects. The fact that we 
are not considering the complete background solution to eq. (|B|) will modify our results only at 
higher orders. 

We first consider the left handed top current J/t = ihl^L- The first non- vanishing term of 
the expansion in eq. (ffl) appears quadratically in the Higgs fields background, %. e. 



i,3 
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m*my), (9) 



*? i =A M =0 



where we take i>£ = {H®, H± , H$, H®*} to be the bubble wall Higgs background. Since we 
are interested in computing the averages in regions where the sphalerons are still active, i.e. 
for Vi/T < 1 [ITJ, higher order contributions to the expansion (0J) can be safely neglected. 



The first non-vanishing contribution to eq. (0) appears only at one loop, and is given by 

W*)> (1) =4 / d A x d*y1m(H°{x)H°{y) m ) G»(x,y,z), (10) 

where h t is the top Yukawa coupling and G^{x, y, z) is the Green function corresponding to the 
diagram in Fig. 1. As we have already remarked, in computing this diagram one must use 
the propagators in the unbroken phase. Note that this contribution vanishes if the phase 82 of 
the Higgs field H® is a constant. The scale of the external momenta p are set by the space- 
time variation of the phase, which, in the case of interest for us of thick bubble walls (L w ~ 
(10 — 100)/T) fll2|, allows us to neglect contributions of order p 2 /T 2 in the high temperature 
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Figure 1: The 1-loop contribution to the left-handed top current. 



expansion. Computing Q )1 (x, y, z) in the MS renormalization scheme at the scale /x, we obtain, 
in the reference of frame of the thermal bath, 

<TO> (1) * tiWMfc) ( lo S [jfj^j + I) > (11) 

where A f = vr 2 exp(3/2 - 2-f E ) ~ 13.944. 

Eq. (|TlD shows the expected dependence on /i 2 and which, in comparison to the 

original result given in ref. 0, gives a suppression factor C(/i t v 2 /7r T) 2 . 

A graph similar to that in Fig. 1 for the right handed top quark leads to a contribution to 
(Jf* (z)) given by (Jj* (z))^ 1 ' = — (J? {z))^ 1 ' . For the other fermion species, one finds analogous 
results, in which ht is replaced by the appropriate Yukawa coupling, and V2(z) (vi(z)) appears 
for the up (down)-type fermions. 

A contribution to (Jf (z)) proportional to Im^H^H^) = v^sinS, and to /i 2 , appears at 
two loops, given by the graph in Fig. 2. Since the computation has to be performed in the 
unbroken phase, we must use resummed propagators for the Higgs fields in order to deal with 
the IR divergences ]13| . 



In the unbroken phase the Higgs spectrum contains two complex neutral fields and two 
charged fields. The resummation can be achieved by considering the propagators for the eigen- 
states of the thermal mass matrix, which have masses given by 

M 2 2(T) = m\(T) + m 2 (T) =F y/jmgj - m 2 (T) 2 ) 2 + 4mf (T) (12) 

where the mf(T) are the thermally corrected mass parameters of the potential (|7|), m 2 (T) ~ 
m 2 + 3g 2 T 2 /16, m|(T) ~ m\ + /i 2 T 2 /4, while m\{T) receives only logarithmic corrections in T, 
which were computed in ref. ||. 



5 




Figure 2: The 2-loop contribution to the left-handed top current. The scalar internal lines can 
be either neutral or charged Higgs fields. 

Assuming p 2 <C M 2 2 (T) <C (27rT) 2 , where p is again the external momentum p ~ 1/L W , it is 
straightforward to see that the non-zero Matsubara modes of the bosonic integral are strongly 
suppressed. Keeping only the zero mode, we obtain 

<*m> w - -^ ju^mw ( log + 5 a " w *- w siaS[z)] ■ (13) 

As for the one- loop result, the two loop contribution with neutral Higgses in the internal lines 
for (Jt R (z)) is opposite to the one for (J^ L (z)). This is no longer true when the graphs with 
charged Higgses in the internal lines are taken into account. In such a case, (Jf* (z)) gets a a 



contribution equal to eq. (]T3[) , whereas the result for {J^ L (z)) is analogous but proportional to 
hi instead of h\ . The charged Higgs loops give also rise to a non- vanishing left handed bottom 
density opposite to eq. fllTf). As we will discuss in the following, this fact may have interesting 
implications for the spontaneous baryogenesis mechanism. 

On the background of the CP violating bubble wall non-vanishing Higgs currents 
Jh. = i (^HjV^Hi — V^HjHij, where T>^ is the covariant derivative, are also present. In the 
case of the neutral Higgses a contribution already appears at the tree level, see Fig. 3, and is 
given by the classical current 

<J£oW> {3) = -2v 2 (z)d^(z). (14) 

i 

Note that, again, the individual Higgs currents are zero if the phases are constant. Moreover, as 
one can easily see from the equation of motion (|8]), the total Higgs hypecharge current vanishes 
at the tree level. 

At one loop, Fig. 4, we have a further contribution to the neutral Higgs current and the 
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Figure 3: The tree level contribution to the neutral Higgs current. 
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Figure 4: The 1-loop contribution to the neutral and charged Higgs currents, 
first non- vanishing one to the charged Higgs current. For the neutral Higgs currents we obtain 



Each charged Higgs gets a contribution equal to that of the neutral Higgs belonging to the 
same doublet. As expected, also the Higgs currents vanish in the limit of vanishing Vi(z). 

In summary, we have made use of a field theoretical approach based on a expansion in the 
background fields around the unbroken phase to compute in a consistent way the perturbations 
to the particle currents induced by a CP violating bubble wall background. We have shown 
that the various contributions arise at the tree level, in the case of the Higgs currents, or as 
loop effects. In this way it has been possible to avoid the various ambiguities inherent to the 
traditional approach based on the rotation of the fields, recovering all the expected suppression 
factors and including the Higgs fields. 

In order to get a feeling of the implications of these results for the spontaneous baryogenesis 
mechanism, we can first consider an adiabatic approximation similar to that discussed in ref. 



<j£oW> (4) = (r H o(z)) 



A 5 TmljT) 
8tt (M 1 (T) + M 2 {T)) 



3^ l Vl (z)v 2 (z) sin 6(z)\, 



(15) 
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II]. Inside the bubble walls the (B + L) -violating sphaleron transitions are biased since a 



non-zero local equilibrium value for (B + L) is induced. (B + L) eq is a linear combination of 
the expectation values of all the currents conserved by the interactions in equilibrium inside the 
bubble wall. Assuming that gauge flavor diagonal, top Yukawa and Higgs-Higgs interactions 
are in equilibrium Q, the conserved charges are Q', (B — L)', BP 1 — B% — 1/2 (Bi + B%) and 
Y' = YH+Y tL +Yt R +Yb L + l/3Y l ' ep where the prime means that only particles in equilibrium must 
be considered, and Y/ ep is the hypercharge of the leptons in equilibrium. The values of these 
conserved currents are made up by the background perturbations computed above. Neglecting 
the contributions to the fermionic currents not proportional to h 2 , it is straightforward to see 
that the contributions to (B + L)eq from the two Higgs currents cancel each other both at 
the tree level and at one loop, and among the conserved charges only Y' gets a non-vanishing 
contribution at one loop. The two loop diagrams where charged Higgses are exchanged give a 
further contribution to Q' and Y' . 

The local equilibrium value (B + L)eq enters a rate equation of the form 
B ~ — V sp{B + L)eq/T 3 where Tsp is the rate of the sphaleronic transition, describing the 
generation of baryon number inside the bubble walls. As a consequence, the force driving 
baryon number violation is just a one loop effect. 

Regarding the final value for the baryon asymmetry in the adiabatic approximation, the 
suppression terms that we have found would give rise to a suppression factor 0(h 2 v 2 /ii 2 T 2 ) 
where v co is the value of the Higgs fields for which the sphaleron transitions cease to be effective. 
Following ref. [jllj and taking v co /T ~ a w /g we see that typically a suppression 0(1CT 4 ) arises. 
In this case, it might be hard to reconcile the observed value for the baryon asymmetry with 
this mechanism for baryogenesis in the adiabatic limit, both in the case of spontaneous and of 
explicit CP violation in the Higgs sector. 

Nevertheless, it has been recently shown that particle diffusion may play an important 
role in the description of the spontaneous baryogenesis mechanism. As a consequence of the 
perturbation induced by the bubble wall background, and of the different diffusion coefficients 
of different particle species, asymmetric densities are formed not only inside the bubble walls, 
as in the adiabatic approximation, but also in front of it, where the sphaleron transitions are 
not suppressed. These asymmetries are then transformed into a baryon asymmetry mainly in 
the region in front of the bubble wall, in a scenario similar to that occurring in the case of thin 
bubble walls. 

In order to improve the adiabatic approximation a system of kinetic equations describing 

2 In this approximation we also assume that the hypercharge violating processes whose rates are suppressed 
by powers of v(T)/T are out of equilibrium. 
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particle interactions and diffusion should then be solved, along the way of, e.g., ref. 0. In this 
case, the quantities that we have computed would represent the sources of the departure from 
thermal equilibrium. This task goes beyond the scope of this letter, however an important, 
and maybe helpful, difference with respect to the situation considered in ref. may be 
pointed out. As we have already discussed, in that paper the perturbation to the i-th particle 
species was simply proportional to y % d^5 T 2 , where y % is the hypercharge of the particle. This 
proportionality was due to the fact that a hypercharge rotation has been made to remove the 
Higgs phase from the Yukawa couplings. Such a dependence does not appear in our results, 
neither for the fermions nor for the Higgses. Instead, as long as the loops with charged Higgses 
are neglected, we obtain that, e.g., (Jt L ) = ~(Jt R )- This behaviour reflects the fact that in the 
Yukawa couplings left and right handed quarks have opposite couplings to the imaginary part 
of the Higgs field. When charged Higgses are taken into account, the contribution to (Jt R ) at 
two loops is opposite to that to (J{? ) while (Jt L ) and (J^ ) get also opposite contributions, but 
only proportional to h 2 . Concerning the Higgs currents, they are proportional to A5 instead 
of h 2 and moreover, from eq. (|T^) we see that they are the same for the two neutral Higgses, 
which have opposite hypercharges. 

In ref. |§ it was noted that, since all the perturbations were proportional to the hypercharge 
(and to 5), the only processes which could be biassed were the hypercharge violating ones, 
which have rates suppressed by powers of Vi(z). Due to the smallness of these rates on the 
outer edge of the bubble wall, this gave rise to an important suppression on the value of the 
fluxes of particles diffused in the unbroken phase. On the other hand, our results do not 
exhibit such proportionality, then also the (unsuppressed) hypercharge conserving processes 
are actually biassed. As a consequence, inserting our perturbation in the kinetic equations one 
might expect that in the unbroken phase fluxes of roughly the same magnitude than those of 
ref. would be found. A numerical treatment of the problem using the results of this letter 
and properly treating baryon number violation in the unbroken phase would then be highly 
desirable in order to understand whether diffusion can really help in overcoming the pessimistic 
results of the adiabatic approximation. 

It is a pleasure to thank W. Buchmiiller, J.R. Espinosa, F. Feruglio, and J. Ignatius for 
useful discussions. 
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